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Abstract 



We present the cross section for pair production of scalar leptoquarks LQ in pp 
collisions, p + p ^ LQ + LQ + X, at the Fermilab Tevatron in next-to-leading order 
5^ \ QCD. Including the higher-order corrections stabilizes the theoretical prediction and 

increases the size of the cross section for renormalization/factorization scales close 
to the mass of the leptoquarks. This leads to an increase of the lower bound on the 
mass of scalar leptoquarks by up to 15 GeV with respect to earlier analyses. 



1. Theoretical speculations have focussed on two different elements to explain the 
recently observed surplus of DIS e^p events at HERA contact interactions at an 
effective scale Aeq ^1-5 TeV, and narrow resonance formation at a mass scale M ~ 
200 GeV. The resonance interpretation is based in particular on the HI data which appear 
to cluster in a narrow range at invariant (eg) masses of about 200 GeV. Such resonances 
1^ can be identified with scalar squarks in supersymmetric theories with i?-Parity breaking 
or with leptoquarks in general. 

Leptoquarks encompass a large variety of particles which can be classified in multiplets 
according to spin, isospin and hypercharge . The production of leptoquarks in positron 
collisions with sea quarks requires large Yukawa couplings, which can in general not be 
reconciled with bounds from rare decay processes, nor with the early e~p scattering data 
from HERA. However, the Yukawa couplings in positron collisions with valence quarks, 
A ~ 1/30, are so small that they are not in conflict with any such bounds. Moreover, they 
cannot be ruled out by present limits on eeqq contact interactions from LEP ^ or from 
the Tevatron [T^. The most powerful competitor in this scenario is the pair production 



of leptoquarks ||TT| at the TevatronQ 

p + p^IQ+lQ + X (1) 
Mass bounds on vector leptoquarks, which decay solely into charged leptons and quarks. 



appear significantly above 200 GeV in the analyses of Refs. |ll2| , p!3| . Even if the unknown 
anomalous couplings of vector leptoquarks are chosen such as to minimize the cross sec- 
tion, the minimum mass is 210 GeV, see MM- The present limits for scalar leptoquarks 



are near 200 GeV |jT2| . The experimental bounds will be refined further by the two Teva- 
tron experiments in the near future. The limits for scalar leptoquarks are particularly 
important since they are essentially parameter- free. The cross sections for pair production 
of scalar leptoquarks involve only the strong coupling constant, and they do not depend 
on unknown Yukawa couplings. 

Anticipating the refinement of the Tevatron limits on the masses of scalar leptoquarks, 
a solid theoretical prediction of the production cross section^ is mandatory, the more the 
gap narrows between the Tevatron mass bounds and the HERA mass estimates. Based 



on previous experience from the production of squark pairs in hadron collisions |T5 
it is expected that higher-order QCD corrections increase the production cross section 
compared to the predictions at the Born level. Experimental mass bounds are therefore 
shifted upwards. Moreover, by reducing the dependence of the cross section on spurious 

^Another source of leptoquarks at the Tevatron is associated production pp — > LQ + 1 + X The 
cross section for this process depends on the leptoquark Yukawa couphng; it is significantly smaller than 
for leptoquark pair production j^. 

^Soft gluon corrections to the production of leptoquark pairs have been discussed in Ref.||l^ [based, 
though, on erroneous Born calculations]. 
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parameters, i.e. the renormalization and factorization scales, the cross sections in next- 
to-leading order QCD are under much better theoretical control than the leading-order 
estimates. The next-to-leading order analysis will be presented in this letter. 

2. The basic processes for the production of leptoquark pairs at the Tevatron are 
quark-antiquark annihilation and gluon-gluon fusion: 

q + q-* LQ + lQ 

g + g^IQ + TQ (2) 

Any non-pointlike structures of leptoquarks which may occur in compositeness scenarios, 
are expected at scales only above ~ 1 TeV; in other scenarios leptoquarks are generically 
pointlike particles. The gluon-leptoquark interactions are therefore determined by the 
non-abelian SU(3)c gauge symmetry of scalar QCD so that the theoretical predictions for 
the pair production of scalar leptoquarks are parameter-free. 

The diagrams corresponding to the processes (Q), are shown in Fig.|I|a; the only new 
element of scalar QCD is the quartic coupling between gluons and leptoquarks which 
follows from the SU(3)c gauge invariance of the interaction. The cross sections of the 
parton processes (EI) may be written as |]TB 



,-.2 



96s 



p (41 - 31/32) + (I8/32 - /^^ - 17) log - 



(3) 



where is the invariant energy of the subprocess and {3 = y 1 — AM^q/s. The cross 
sections coincide with the cross sections for squark-pair production in the limit of large 
gluino masses p!5| , p!7| . The quark-antiquark annihilation is the driving mechanism at the 
Tevatron for large leptoquark masses. 

The QCD radiative corrections to the order ctg include virtual corrections, the brems- 
strahlung of gluons and contributions from gluon-quark collisions. The virtual corrections 
can be classified in self-energy diagrams and vertex corrections for quarks, gluons and 
leptoquarks according to the rules of standard and scalar QCD, and initial/final state 
interactions. In addition to gluon radiation off all colored lines, gluons can also be emitted 
from scalar vertices. Finally, the inelastic Compton process, diagram (^), must be added 
in order of the cross section. The amplitudes have been evaluated in the Feynman 
gauge. After the singularities are isolated by means of dimensional regularization, the 
renormalization has been carried out in the MS scheme. The masses of the light quarks 
(u,d,s,c,b) have been neglected while the mass parameter of the leptoquark has been 
defined on-shell. We have chosen a renormalization and factorization scheme in which the 
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massive particles (top quark, leptoquark) are decoupled smoothly for momenta smaller 
than their mass |T^. This implies that the heavy particles do not contribute to the 
evolution of the QCD coupling and the parton densities. The computation of the cross 
section for gluon emission has been performed by adopting the phase space slicing method 
(see e.g. fl^): a cut-off A has been introduced for the invariant mass of the leptoquark- 
gluon system in the final state, which separates soft from hard gluon radiation. If both 
contributions are added up, any A dependence disappears from the total cross section for 
A — > 0. The infrared singularities cancel when the emission of soft gluons is added to the 
virtual corrections. The remaining coUinear initial-state singularities are absorbed into 
the renormalization of the parton densities |^D[, defined in the MS factorization scheme. 

The perturbative expansion of the total parton cross section can be expressed in terms 
of scaling functions. 



(4) 



with z, j = g,q,q denoting the initial-state partons. For simplicity, we have identified the 
renormalization scale with the factorization scale fiR = fip = fJ'- The scaling functions 
depend on the invariant parton energy v'S through rj = s/AM^q — 1 and, very mildly, 

on the ratio of the particle masses = rritop/MiQ. The scaling functions f-^'^^ and f^^^^ 
are displayed in Fig.|| for the quark-antiquark, gluon-gluon and gluon-quark channels, 
respectively. The scaling functions f^p are decomposed into a "virtual + soft" (V+S) 
part, and a "hard" (H) gluon-radiation part; the In-' A (j = 1, 2) singularities of the (V+S) 
cross section are mapped into (H), cancelling the logarithms so that these functions are 
independent of A in the limit A ^ 0. 

From Fig.|^(a) we can infer that the next-to-leading order corrections to the gluon- 
gluon channel are very important near the threshold \/§ > 2Miq. The large size of the 
corrections close to the threshold is partly due to gluon radiation off the initial-state par- 
tons, which generates contributions of the type pin-'P (j = 1,2); the In^P terms are uni- 
versal and can be exponentiated [^. At the threshold, 2Miq, the next-to-leading 
order cross section for the gg initial states is non-zero: The Sommerfeld rescattering con- 
tribution, due to the exchange of Coulomb gluons between the leptoquark pair in the 
final state, gives rise to a 1//3 singularity which compensates the phase space factor 
At large parton energies, the hard coefficients /gg and fg]^ approach plateaus, that are 
built up by the flavor-excitation and gluon-splitting mechanisms. The exchange of gluons 
in the t- and u- channels leads to an asymptotically constant cross section, which is to be 
contrasted with the scaling behavior ~ 1/s of the leading-order process. 

The pp cross section is found by folding the parton cross sections with the gluon and 



4 



light-quark luminosities in pp collisions: 

dT^aij{s = Ts) (5) 

Due to the dominating qq luminosity for large parton momenta, the total cross section is 
built up primarily by quark-antiquark initial states for leptoquark masses Miq ^ 100 GeV. 
For the numerical analysis we adopt the CTEQ4M parametrization of the parton densities 
p^ . The QCD coupling is evaluated in the MS scheme for nif = 5 active flavors and 



A^^) = 202 MeV, and the top quark mass is set to mtop = 175 GeV |p3 |. 

The scale dependence of the theoretical prediction is reduced significantly when higher 
order QCD corrections are included. This is demonstrated in Fig.^ where we compare the 
renormalization/factorization scale dependence at the leading and next-to- leading order 
of the total cross section. For a consistent comparison of the LO and NLO results, we 
have calculated all quantities [i.e. as{fi'^), the parton densities, and the partonic cross 
sections] in leading and next-to-leading order, respectively. The scale dependence of the 
leading-order cross section is steep and monotonic: Changing the scale from n = 2Miq to 
/i = Miq/2, the LO cross section increases by 100%. At next-to-leading order the scale 
dependence is strongly reduced, to about 30% in this interval. The NLO cross section 
runs through a broad maximum near n ~ MiQ/2, which supports the stable behavior 
in /i. 

3. The QCD radiative corrections enhance the cross section for the production of 
leptoquarks above the central value ^ ~ ^iq- If? as often done, a LO cross section 
is calculated in a hybrid form, i.e. using parton cross sections in Born approximation 
but NLO parton densities and as in two-loop order, the enhancement is close to 40%. 
If, moreover, the so-defined LO cross section is evaluated at large scales ~ V^, the 
enhancement in next-to- leading order is as big as ~ 70%, nearly independent of the 
leptoquark mass. The convergence of the perturbative approach should however be judged 
by examining a properly defined i^'-factor, K = cttsslo / (^lo , with all quantities in the 
numerator and denominator calculated consistently in NLO and LO, respectively. In 
the interesting mass range between 150 < Miq < 250 GeV, these K-factors vary only 
between 1.20 and 1.08, as shown in Table |l]. They are small enough to assure a reliable 
perturbative expansion.^ 

■^The results nearly coincide with the cross sections for the production of squark-antisquark pairs in 
the limit of large gluino mass. Supersymmetry predicts quartic self-interactions to order between the 
squarks. No such self-couplings have been considered hitherto for leptoquarks in general. These couplings 
affect the cross sections only through rescattering corrections involving heavy leptoquark loops so that 
their impact on the cross section is small. We have checked that the results after removing these diagrams 
are identical for large gluino masses. 
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The impact of the next-to-leading order QCD corrections on the present experimental 
lower mass limits for scalar leptoquarks is illustrated in Fig.^. We compare the NLO 
result based on the default settings (CTEQ4M parton densities and fi = Miq) with the 
LO cross section adopted in earlier analyses (CTEQ3M parton densities and two-loop 
cts, evaluated at /i equal to the partonic centre-of-mass energy \fl ).0 Taken at face value, 
the next-to-leading order corrections increase the mass limit for a first-generation scalar 
leptoquark by about 15 GeV. [The shift is much smaller if the LO cross sections are 
evaluated for LO parton densities and as at the renormalization/factorization scale /x = 
MiQ, as demonstrated by the broken line in Fig.^.] The shaded band reflects the remaining 
theoretical uncertainty at NLO due to the choice of the renormalization/factorization 
scale when is varied in the range Miq/2 < < 2Miq. Since the cross section in 
the interesting mass region Miq ^ 150 GeV is built up mainly by the quark-antiquark 
channels, thus based on well-measured parton densities, the variation between different 



parton parametrizations (CTEQ4M [||, GRV |26| and MRS(R2) ^) is less than 5 %. 

Evaluated on the basis of the NLO cross sections presented in this letter, the data 
from the Fermilab Tevatron, which were presented in Refs.|jl2|, lead to parameter-free 
lower limits of about Miq > 190 GeV and 210 GeV for scalar leptoquarks decaying to 
charged leptons in the present DO and CDF analyses. 
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MiQ [GeV] 


(Jqq CFgg Ggq (Jtot [pb] 


K 


150 


LO 0.741 0.244 0.985 
NLO 0.722 0.490 -0.028 1.184 


1.20 


175 


LO 0.318 0.071 0.389 
NLO 0.311 0.146 -0.010 0.447 


1.15 


200 


LO 0.142 0.022 0.164 
NLO 0.141 0.047 -0.004 0.184 


1.12 


250 


LO 0.030 0.003 0.033 
NLO 0.030 0.006 -0.001 0.035 


1.08 



Table 1: LO and NLO results for the total cross section p + p ^ LQ -\- IQ + X at the 
Tevatron energy y/s — 1.8 TeV for various values of the leptoquark mass Miq. All quan- 
tities [as{Q'^), the parton densities, and the partonic cross sections] have been calculated 

consistently in leading and next-to-leading order, respectively. Also shown is the K factor 
defined as K — (Jj^io/ cfiq- The CTEQ4M(L) parton densities with the associated values 
of and the central scale fi = Miq have been adopted. [The negative sign of dgq is a 
mere artifact of subtracting collinear initial-state singularities via mass factorization.] 
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Figure 1: Generic diagrams for pair production of scalar leptoquarks in hadron collisions: 
(a) qq annihilation and gg fusion; (b) the gluon-quark subprocess. 
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Figure 2: The scaling functions for scalar leptoquark pair production in gg (a), qq (b), 
and gq (c) collisions versus 7] = s/AMIq — 1. The notation follows Eq.^; B denotes 
the lowest-order scaling function, V+S the sum of virtual and soft corrections, H the 
contribution of hard gluon emission, SC the scale- dependent scaling function. 
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Figure 3: Renormalization/ factorization scale dependence of the total cross section p+p 
IQ + IQ + X at the Tevatron energy a/s = 1.8 TeV. Parameters as described in the 
text. The arrow indicates the value of the average invariant energy <s>^/^ in the hard 
subprocess. 
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Figure 4: The total cross section for the production of leptoquark pairs, p + p ^ IQ + 
IQ + X , at the Tevatron energy y/s = 1.8 TeV as a function of the leptoquark mass Miq. 
The next-to-leading order result (NLO) is compared with various options of parameters in, 
partially hybrid, LO calculations. The variation of the NLO cross section with the value 
of the renormalization/ factorization scale is indicated by the shaded band. The increase of 
the lower bound on the leptoquark mass, compared to the earlier analyses, is demonstrated 
by the horizontal arrow. 
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